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Abstract 

Q \ We study the four-body exclusive semileptonic baryonic B decays of B~ — )• ppi~U£ {£ = e,/i, r) 

in the standard model. We find that their decay branching ratios are about (1.0, 1.0,0.5) x 10~^, 
respectively. In particular, the electron mode is close to the corresponding CLEO's upper limit 
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r~| . of 5.2 X 10 ^, while all results are about one or two orders of magnitude larger than the previous 
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estimated values for the inclusive modes of i? — > BB'^i/. Clearly, both B-factories of Belle and 
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BaBar should be able to observe these exclusive four-body modes. 



I. INTRODUCTION 

In the semileptonic B — )■ Miij£ decay with a meson M and a charged lepton i, the iue pair 
involves no direct QCD interaction so that the theoretical description of the amplitude can be 
reduced to a simple form with the 5 — )■ M transition. For example, the rate for B^ — )■ ir'^e'h'e 
is proportional to |Vub/+(g^)p, where the form factor f+iff) foi' the B^ — )■ 7r+ transition 
depends on the momentum transfer squared, q^ . This benefits the precision measurement 
of \Vub\i where \Vub\ is one of the least known Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements [l|, |2| in the Standard Model (SM). As long as we choose a point q^ = qf in the 
decay spectrum, the corresponding data point with other parameters can be fixed to extract 
the value of \Vyh\. However, f+io^) relies on the calculations in the QCD models, such as 



quark models 



3|, lattice QCD j^, and Light Cone Sum Rules 



5|. Starting with qf and 



Vub\, one is allowed to inversely extract the g^ dependence of f+{q^) in different q^ intervals 



from the measured data 



-[id]- The extraction compared with various theoretical models 



hence improves the knowledge of /+(g^). Moreover, such extraction also provides crosschecks 
for the i? — 7- p and B — )■ rf^'' transition form factors 



-Il0|. In particular, the size of the 
gluonic singlet contribution llNl3| to the B ^ rj' transition to explain the unexpectedly 



large two-body hadronic B — )■ Kr]' decay rate has been constrained by measuring B — )■ rf^'Hy 



decays 



10|. Similar to the mesonic cases, it should be interesting to extend the study to 



baryonic decay modes, such as i? — )> BB'£z/ with BB' being a baryon pair, to investigate 
the B — )■ BB' transition form factors, which have been used as theoretical inputs in the 
three-body B — t- ppM decays. 

The factorizable amplitudes for the three-body baryonic B — )■ BB'M decays are normally 
classified into current and transition parts, given by 



Acx(BB'|(gig2)|0)(M|(g36)|S) 
^roc(M|(gig2)|0)(BB'|(g36)|fi) 



(1) 



respectively, where {qiq2) and {q^h) stand for the weak currents. The matrix elements of 
— )■ BB' in Ac are presented as the timelike baryonic form factors, for which the theoret- 
ical calculations are available, such as the approach of the pQCD counting rules 14l-ll6|. 
Consequently, the observed branching ratios for 5° — )■ npD*^ 17|, B~ — )■ Apvr" 



B -^ AAA'(*) 



22 



2lj and 
23| can be explained due to the v4c-like amplitudes |24J-|32| . On the other 
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FIG. 1. Baryonic B decays with the B ^ pp transition, where (a) B 
(c) B -^ ppK^*\ and (d) B^ -^ ppD'^*'>°. 



pplvf,, (b) B — )• ppiT 



land, the measured decays of 5 — )> ppK^*\ B~ 
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391], shown in Fig. [H are considered to have Aj- as their amphtudes 25| 
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42| . To explain the data, the transition matrix elements of i? — )■ BB' are parameterized in 
terms of various form factors 
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42| . As seen in Fig. [H the decay of -B — > ppiue 



is also a ^7- type as -B — ?■ ppM. It is clear that the observation of B^ — )■ ppiug shall directly 
determine the transition form factors, which have been widely used to explain the B — ?■ ppM 
data as theoretical inputs. In analogy with the timelike baryonic form factors, similar mo- 
mentum dependences of the transition form factors may be chosen, which can be justified 
by investigating the shape of the invariant mass spectrum for the B~ — ?> ppi'Ui decay. 
Moreover, we expect that the measurements of the angular distributions in B^ — )■ ppe~Ue 



will provide some information to understand the unexpectedl 
asymmetries of A0{B~ — )• ppK~) = 0.45 ~ —A0{B~ -^ ppix") [35|. 



arge angular distribution 



At present, the CLEO Collaboration has given an experimental upper limit: 43|] 

B{B~ -^ ppe-Ue) < 5.2 x 10~^ (90% C.L.) , (2) 

BB'iu decays 



while the theoretical estimation has only been done for the inclusive B 
with charmless dibaryons, given by 4^ 

-5 



B{B -^ BB'£z/) 



10' 



10' 



(3) 



In this paper, we concentrate on the exclusive four-body semileptonic baryonic decay of 
B~ — )■ ppiue {i = e, fi, or r). In particular, we will study its decay branching ratio in the 
SM. 

The paper is organized as follows. In Sec. II, we provide the formalism, in which we show 
the decay amplitude and rate of B" — )■ ppiu£ along with the definitions of the transition 
form factors of i? — t- BB'. We give our numerical results and discussions in Sec. III. In Sec. 
IV, we present the conclusions. 



II. FORMALISM 



In terms of the effective Hamiltonian, given by 



U{h ^ uiu) = ^^ n7^(l - 75)6 £V(1 - 75)^ . 



(4) 



for the h ^ u transition with the W boson emission to dv at the quark level, we easily 
factorize the amplitude for the B~ — > ppCun decay to be 



A{B- -> pplv,) = ^^{pp\u^^{l - ^,)h\B-) £V(1 - ^,)v, 



(5) 



where we have parameterized the amplitude in terms of the transition matrix element of 
B — )■ pp. With Lorentz invariance, the most general forms of the B — )■ BB' transition form 
factors can be written as 31 1 



(BB'|g'7^6|5) = m(pB)[fi'i7M + 92^(r^^uP'' + QsPf, + 9a{p^' + Pb)m + 9b{Pw - Pb)/.]75^'(Pb')> 
(BB'|g'7^756|S) = w(pb)[/i7m + /2^^m^p'' + /aP^ + /4(Pb' + Pb)m + /sIPb' - Pb)m]^(Pb').(6) 

with p = Pb — Pb — Pb' fo^^ the vector and axial- vector quark currents, respectively. For the 
momentum dependences of fi and Qi, we can rely on the results in the B -> ppM decays 
as they share the same B — )■ BB' transition form factors. Since the pp invariant mass 



distributions for B — t- ppM have been observed to peak near the thresho 



out at the large energy region, inspired by the pQCD counting rules |l4|-|l6l. l27l |. we simply 



d area and flatten 



take the form factors as 



4Q| 






(7) 



with n = 3 and t = [pp + ppY = i^pp where Df. and Dg- are constants determined by the 
B — )■ ppM data. Note that the number of n = 3 is for three hard gluons as the propagators 
to form a baryon pair in the approach of the pQCD counting rules, where two of them attach 
to valence quarks in pp, while the third one kicks and speeds up the spectator quark in B. 
In terms of Eqs. ([5]), ([6]), and ([7j), the amplitude squared |^p by summing over all fermion 
spins becomes available. 

We then need the kinematics for the four-body B{p§) -^ B(pb)B'(pb/)£(p^)i/(Pp) decay 
to integrate over the phase space. As the formalisms in K14, D14, and B14 decays given in 



Refs. 



45l447l|. we use five kinematic variables, s = {pi + Pf^ 



ra 



ipy 



t, ^B, 6'l, and to 



describe the decay. The three angles 9b, Ol, and are depicted in Fig. [21 where the angle 



^B(L) is between pb (p^) in the BB' (iu) rest frame and the hne of flight of the BB' (£p) 
system in the rest frame of the B meson, while the angle (p is from the BB' plane defined 
by the momenta of the BB' pair to the iu plane defined by the momenta of the iu pair in 
the rest frame of B. The partial decay width reads 




FIG. 2. Three angles of Ob, Oi,, and (p in B ^ BB'ii?. 
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where X, /3b, and /3l are given by 



X 



-{ml-s-tY-st 



1/2 



/3l = -A^/^(s,m^,mp), 



(8) 



(9) 



respectively, with A (a, b, c) = a^ + b"^ + c^ — 2ab — 2bc — 2ca. The regions for the five variables 
of the phase space are given by 



"u) ) 



(m£ + mp)^ < s < (m^ — Vt)^ , {rriB + m-B')^ < ^ < (m^ — m^ — rrii 

0<^L, ^B <7r, O<0<27r. (10) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

In our numerical analysis, we take \Vuh\ = (3.89 ± 0.44) x 10"'^ from the PDG [48']. To 
deal with Dg. and Df^ in Eq. (I7j), it is helpful to use the approach of the pQCD counting 
rules again, where with SU{3) flavor and SU{2) spin symmetries the vector and axial- vector 



currents are incorporated as two chiral currents in the large t limit [27 



31 



40|. Consequently, 



Dg. and Dj^ from the vector currents are related by another set of constants D\\ and D 



the chiral currents. Exphcitly, for the B -^ pp transition form factors we have 3l|, |40 | 
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(11) 



3 M 3 ir J^ 3 " 3 II' «^ 3 

with j = 2, 3, .., 5, where their values are determined by fitting the data of the total branch- 
ing ratios, invariant mass spectra, and angular distributions measured in the B — )■ ppM 
decays. To illustrate our results, we adopt the values in Ref. |3l| : 

(L)||, Dj^) = (67.7 ±16.3, -280.0 ± 35.9) GeV^ 

(Dl Dl Dj, Dl) = 

(-187.3 ± 26.6, -840.1 ± 132.1, -10.1 ± 10.8, -157.0 ± 27.1) GeV^ . (12) 

Thus, the total branching ratios of B~ — )■ ppiug are found to be 

B{B~ -^ ppe~iJe) = (1.04 ± 0.26 ± 0.12) x 10"^ , 
B{B~ -^ PPli'y^,) = (1.04 ± 0.24 ± 0.12) x 10"^ , 
B{B- -^ ppT-Ur) = (0.46 ± 0.10 ± 0.05) x 10"^ , (13) 



where the two errors in Eq. flT^ are from those in Eq. flT2]) and \Vub\, respectively. The 
invariant mass spectra and angular distributions for i?~ — )■ ppe~Ue are shown in Fig. |3l The 
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FIG. 3. Invariant mass spectra as functions of the invariant masses rupp and rueue and angular 
distributions as functions oi 9i {i = B, L) for S^ — )• ppe~v'e, respectively. 



integrated angular distribution asymmetries, defined by 



Ae 



/o H^dcosO^ - Ti 



dB 



dcosB; 



■dcosOj 



/o l^dcose. + /-I 



° "^ dcose, ' 



B,L) 



(14) 



dcosdi 



are obtained to be 



Ae^iB- -^ ppe^Ue) = 0.06 ± 0.02 , 

Ae^iB- ^ppe-Ue) = 0-59 ±0.02, (15) 



where the errors are from those in Eq. flT2|) . 

Since our result on B{B~ — )■ ppe~Ue) in Eq. flT3|) is around 1.0 x 10~^, which is the same 
order of magnitude as those of the well measured mesonic B decays at Belle and BaBar, such 
as B{B^ — )• 7r^(p'^)i~h'ii) and B{B^ — )• p^i'ue), this four-body mode should be observed at 



these B- factories 49|]. Moreover, as seen from Fig. [3|i, the B^ — )■ ppe^Ve decay inherits the 
same threshold enhancement as those in the three-body baryonic B decays, resulting from 
the adoption of 1/t^ for the momentum dependence in the B^ — )■ pp transition form factors. 
The spectrum in Fig. [3]3 reflects the fact that in the helicity structure the amplitude of the 
e'Ve pair is proportional to (-Eg + E^^). 

It is interesting to note that our study of B^ — )■ pplvf, is similar to that of B^ — )■ 



ppK*~ 40l-l42| . The terms related to g2 and /2 in the B~ — )■ pp transition form factors give 
the main contributions to B~ — )■ ppi~vi. Since the pair of the left-handed electron and the 
right-handed anti-neutrino in the helicity structure behaves as one of the polarization vector 
S-{p) with p = Pi + Pue, leading to e ■ p = 0, the contributions from f^ and g^ disappear. 
Those from f^ and g^ are effectively small due to the tiny \Dn\ ~ 10 GeV^. As the branching 
ratio receives the most contribution near the threshold area, the (75(/5)-accompanied term 
{pp ~ Pp) = {Ep — Ep, Pp — Pp) — )• (0, 0) is suppressed. Moreover, since the terms of g2 
and /2 contain a^^p^ , we have the relation \Dg^i^f^-)p\ ~ 300 \p\ GeV^ > |-D/J ~ 200 GeV^ 3> 
\Dg-^\ ^ 20 GeV^, which explains why (72 and f-z prevail over other terms in the B~ -^ ppii^e 
decay. 

Finally, it is interesting to point out that the angular distribution asymmetries in Eq. f lT^ 
in the B~ — !■ ppiug decay are sensitive to new physics, such as the currents of (V^ + A) and 
{S ± P) beyond the SM. Note that Ae^ = 0.59 in Eq. ([15]) (see also Fig. [3t) is caused by 
the £~z/£ pair of {V — A) in the SM, which forms a polarization vector et{p) to couple to the 
left-handed helicity state of the virtual weak boson W*~ . Therefore, a new physics with the 
(y + A) current, which lets the iug pair to be another polarization state £+(p ), must result 



in the deviation of Aq^ in Eq. (fT5|) . Apart from a direct CP violation 4l||, B — ?> ppli'i 



can easily create T-odd triple product correlations (TPC's) to test direct T violation effects. 



Since the three-momenta of pp and those of iui are not in the same plane, pi ■ Opn x pp) 



can be a nonzero TPC observable. Like the case of sa ■ (pa x Pp) in B^ -^ Ap7r~ 50| with 
sa denoting the A spin, there are other TPC observables sg ■ {pp x pp) and Sp ■ {pp x pp) in 
B^ — 7- pplug. These rich TPC observables are expected to be useful to test new physics in 
the advantage of B of order 10"^ much larger than the sensitivity of 10~^ in the B factories. 



IV. CONCLUSIONS 

We have examined the four-body semileptonic baryonic B decay of B^ — )■ ppiue in the 
SM, which proceeds via b — )■ uiu£ at the quark level. The transition form factors of B~ -^ pp, 
which are well studied in the three-body baryonic B — > ppM decays, play the key role in the 
theoretical calculation. We have found that B{B^ — )■ ppiue) = (1.04, 1.04, 0.46) x 10~^ for 
£ = e, fi, T, respectively, which are just a little below the CLEO's upper limit of 5.2 x 10~^ 
for B{B~ — )■ ppe~Ue) but much larger than the previous estimations of 10^^ — 10^^ for the 
inclusive modes of i? — )■ BB'£z/. It is clear that the four-body decays of B~ — )■ ppiue, in 
particular the light charged leoton modes, should be observed by the B-factories of Belle 
and BaBar as well as future B-factories, such as Super-Belle and LHCb. 
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